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CALIBRATION OF COLOR REPRODUCTION APPARATUS 



FIELD OF THE INVENTION 

This invention relates to color calibration for digital color 
reproduction apparatus such as, for example, a color fax or a printer. 

BACKGROUND OF THE INVENTION 

Color reproduction apparatus renders an image from color signals 
representing a color image description. The input to such a color reproduction 
apparatus is an array of color signals (sets of input code values) communicating 
the desired colors of a corresponding array of picture elements in the image to be 
rendered. 

The reproduction apparatus causes cyan, magenta and yellow, and 
sometimes black, dyes to be laid down on a receiving medium in response to color 
signals (its inherent drive code values) applied to its dye marking or producing 
devices. 

A goal is to make the renderings match an original intent as nearly 
as possible. In order to do so, a transform is employed which maps the received 
input code values, normally in RGB space, to the reproduction apparatus' inherent 
drive code values, so that the system faithfully reproduces the color image 
original. Often the transform is implemented by employing a set of three or four 
one-dimensional calibration lookup tables. The process of deriving such a 
transform is called system color calibration. 

SucKapparatus requires re-calibration from time to time due to 
changes in the apparatuKbehavior, as might be affected by the state of processor 
chemicals if any are used, cohar marking system drifts, and changes in receiver 
characteristics. Differences in receiver characteristics might be caused by 
different recipes, variations in receiversNrom one coating event to another, and 
different surface finishes (i.e., glossy versesWtte finishes) as are commonly 
encountered using photographic receivers. \> 
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Calibration of a color reproduction apparatus begins by specifying 
color aim curves. A color aim curves are a set of three or four separate curves, 
one for each of the red, green blue and possibly black channels and each specifies 
the desired density output as a function of the corresponding input code value. It is 
intended that the three (or four) separate curves are to be used in concert and is 
designed to sample the entire density range. It is commonplace that when the input 
code values for red, green and blue channels are equal, that a visually neutral gray 
is to be rendered. The color aim curve defines both the definition of neutral in 
terms of the relationship between measured densities and the relation between 
(common) input code values and output densities. The entries of the table define a 
relationship of density as a function of input code value, or: 

Density = f(iCV) 

It is further commonplace to develop three such functions for a 3 -color printer, 
one for each color. Thus, the family of functions can be expressed as: 

Density Red = ^(iCV r ), 
Density Green = f^JiCVg), and 

Density BIue = f blue (iCV b ). 
These equations are used to find the correct drive code values for each entry in the 
three calibration lookup tables. 

Often, aim curves are linear. A truly neutral density (gray) 
calibration target (FIG. 1) contains a series of aim neutral (gray) density patches, 
which are encoded as grays ranging from white to black. The input code values in 
the target's description are passed through existing calibration lookup tables to re- 
map them into drive code values that the reproduction apparatus will actually use 
to make a print. If the reproduction apparatus is in calibration, the resultant print 
will have the correct aim densities. Mathematical transforms may be employed in 
lieu of lookup tables. The resulting densities of the rendered patches are measured 
using an optical instrument such as a densitometer; and compared to the aim 
densities. If any of the patch density values are different than aim, being either 
too light or dark, or possess a non-neutral hue, then new sets of calibration lookup 
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table values or transforms are derived using a mathematical method such as 
regression so as to increase or decrease the density of the corresponding color. 

The limitation in this process is that the red density is not just a 
function of the cyan drive code value. Rather, it is a function of the drive code 
values of all the other colors. Albeit, the red density has a strong dependency on 
the cyan drive code value and a weaker dependency on the other drive code 
values. This being the case, if one adjusts, say, the magenta drive code value to 
make a correction to the green density, then the red density will change even 
though the cyan code values do not change. Similarly, changes to the cyan code 
values will produce changes to the green and blue densities. Although this known 
calibration process usually works because the cross-dependencies are small, it is 
often necessary to conduct a series of successive calibrations, each using the 
previous output lookup tables, in order to achieve a final calibration within 
specification. 

DISCLOSURE OF THE INVENTION 

It is an object of the present invention to provide for calibration of 
color reproduction apparatus by analyzing the densities read off of a calibration 
target rendered on the reproduction apparatus and by computing a new calibration 
lookup table for bringing the color reproduction apparatus into calibration such 
that the desired behavior of the reproduction apparatus when it is in calibration is 
described by an aim curve which specifies the desired density output as a function 
of drive code values. 

It is another object of the present invention to provide for 
calibration of color reproduction apparatus by analyzing the densities read off of a 
calibration target rendered on the reproduction apparatus and by computing a new 
calibration lookup table or equivalent thereof for bringing the color reproduction 
apparatus into calibration such that the desired behavior of the reproduction 
apparatus when it is in calibration is described by an aim curve which specifies the 
desired density output as a function of drive code values in a single pass. 
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^ Th^present invention involves the rendering of a target containing 

not only neutral densuyspatches but also a sampling of intentional off-neutral 
patches, produced by selectihgsets of red, green, and blue (and sometimes black) 
input code values for those patchefcs^hich deviate from sets of values that would 
5 produce neutral densities, and then, by analyzing how the measured densities of all 
the patches change with changes in code value?J^<teteimming the optimum set or 
sets of the three (or four) input code values that when^^sed together produce a 
desired set of aim neutral density patches in a single iteratii 

According to a feature of the present invention, calibration is 
10 effected by using the color reproduction apparatus to produce a calibration target 
having a plurality of color patches by using input code values which correspond to 
a sampling of color densities. One subset of the color patches is intended to have 
neutral color density values and another subset of the color patches is intended to 
have non-neutral color density values deviating in their mix of red, green, and blue 
15 from each other. The color densities of the color patches produced on the 

calibration target are measured. Error signals characteristic of differences between 
measured color densities and intended color densities are produced. Finally, the 
converting process of said reproduction apparatus is adjusted as a function of the 
error signals so as to provide output color densities closer to the intended color 
20 densities. 

The invention, and its objects and advantages, will become more 
apparent in the detailed description of the preferred embodiments presented below. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the detailed description of the preferred embodiments of the 
25 invention presented below, reference is made to the accompanying drawings, in 
which: 

FIG. 1 is a schematic diagram showing a near-neutral test pattern 
employed in the prior art; 

FIG. 2 is a schematic diagram illustrating a color reproduction 
30 system according to the present invention; 
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FIG. 3 is a plot of measured density data from a target verses input 
code value plotted with the desired behavior; 

FIG. 4 is a chart of patch code values; and 
FIG. 5 is another chart of patch code values. 

DETAILED DESCRIPTION OF THE INVENTION 

The present description will be directed in particular to elements 
forming part of, or cooperating more directly with, apparatus in accordance with 
the present invention. It is to be understood that elements not specifically shown 
or described may take various forms well known to those skilled in the art. 

Inferring to FIG. 2, a color image, represented in digital form as a 
plurality of image intout code values iCV R , iCV G , iCV B , corresponding 
respectively to desired tensities D R , D G , and D B , is inputted by a digital image 
processor 12. The image processor or the reproduction apparatus contains three 
separate calibration lookup tables (or channels) in a 3 -color reproduction 
apparatus (one lookup table foneach color red, green and blue) for converting 
received digital input code values 1CV R , iCV G , iCV B into uniquely-associated 
drive code values dCV R , dCV G , and oCV B , respectively, for a reproduction 
apparatus 14. The image processor or tne reproduction apparatus may contain 
four lookup tables or channels for a 4-colorsreproduction apparatus (black and the 
other three colors). Mathematical transforms n^ay be employed in lieu of lookup 
tables. Reproduction apparatus 14 uses the code values to produce color 
separations for the final rendering. \^ 

Calibration of the color reproduction apparatus begins by 
specifying one or more&im curves of the desired density output as a function of 
input code values for the entirb^iensity range. In FIG 3, the diagonal line 
represents the desired behavior of mb^ystem for one possible aim curve 
description whereby input code values areHinearly related to desired output 
densities. A calibration target containing a seriefes^f both neutral and non-neutral 
aim color density patches is digitized into a series of oeqsity signals and passed 
through image processor 12 to remap the image's input codfe^^alue signals into 



drive oqde values used by reproduction apparatus 14 to make a rendering. The 
resulting obesities of the rendered patches are measured using an optical 
instrument sucn^a densitometer 18; and compared to the intended color 
densities. In FIG. 3, the data points are the actual measured densities from the 
target plotted against the target input code values. If, as in FIG. 3, any of the 
patch density values are differerHsttian intended, being either too light or dark, or 
possess a non-neutral hue, then, usin^a mathematical method such as regression, 
new sets of lookup table values or transfer^ are derived so as to increase or 
decrease the density of the corresponding colorS 

The present invention, recognizes that in the calibration processes, 
each color density is not just a function of the drive code value for that color, but 
rather, it is a function of the drive code values of all the other colors. For 
example, color reproduction apparatus adapted to expose color silver halide paper 
with red, green, and blue light will experience dye sensitivity crosstalk wherein 
the behavior of each color emulsion layer of the receiver is not independent of the 
exposure of the other layers. That is, green light will not only expose the green 
sensitive emulsion layer, but will have some effect on the red and blue sensitive 
emulsion layers as there are overlaps in sensitivity. Accordingly, changing the 
drive code values for one color will result in small effective changes in the other 
color channels. 

Rather than rendering only patches having neutral code values, the 
present invention provides for producing a series of additional patches that vary 
the red, green and blue input code values somewhat from neutral values. For 
example, one might render one patch with neutral input code values, another patch 
in which the red input code value is increased slightly, still another patch in which 
the green input code value is increased, and yet another patch in which the blue 
input code value is increased. In this way, one not only can determine if, and by 
how much the neutral patches differ from neutral, but also how much the densities 
change for known changes in the individual drive code values which are 
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determined by passing the known input code values through the calibration look 
up tables. 

*ere are a number of mathematical paradigms for using this 
information. In a veiysimple example, a certain neutral aim color should be 
achieved with drive code values of, say, 100, 100, and 100. If one renders three 
patches shown in FIG. 4, one ^ith drive code values 104, 99, and 99; the second 
with drive code values 99, 104, 99; ahd the third with drive code values 99, 99, 
and 104, and if in the domains 99 to 104 the behavior is assumed to be linear, 
relations between drive code values and resultihg densities can be written as: 

D r = a*dCV r + b*dC^+ c*dCV b , 
D g = e*dCV r + PdCV g + ^dCV b , and 
D r = i*dCV r 4- j *dCV g + k*d( 
ice there are three sets of drive code values and three sets of 
resulting density pahrts, (D r , D g , D r ), there are nine equations and nine unknowns. 
This is a classical MatmNAlgebra problem with known, published solutions. 
Once the nine coefficients arestetermined, the equations can then be used to solve 
for the unique set of dCV r dC V g aCV b of drive code values that produce the set of 
aim neutral densities, D r D g D b . Givent&e information set forth above, numerous 
schemes using different numbers of variatidft^atches and improved accuracy can 
be thought of by those skilled in the art. The mfhqnum number of patches is 
three. 

would be two problems with the very simple example 
described above. One'problem is that, depending on how the three sets of patches 
are chosen, one may never produce a neutral patch so one never can be sure that 
the neutral is in fact achieved. Th§ second problem is that the functional behavior 
must be not only linear, but the fimctibn must pass through the origin. A simple 
solution to both problems in the above scheme is to add the neutral patch (100, 
100, 100) as shown in FIG. 5, and to add an offsfct to the equations, such that: 



D r = C r + a*dCV r + b*dCV g + c J HCV 1 



D g = C g + e*dCV r + PdCV g + g*dC W and 



N^Db = C b + i*dCV r + j *dCV g + k*dCv b . 

Adding yet another p^bi^total of 5 patches) would allow a quadratic term in the 

most sensitive channel. Equatkms might look like: 

D r = C r + a*dCV r +^^CV g + c*dCV b + p*dCV r 2 , 
D g = C g + e*dCV r + PdC\N^g*dCV b + p g *dCV g 2 , 

D b = c b + i*dcv r + j*dcv g + k^iqv b + p b *cv b 2 . 

Adding additional patches would allow additional degfcs^s of freedom to the 
equations and/or regression analysis to reduce the effect ofnqise. 

\ The collection, called "cluster," of patches including the neutral 
patch and the deviated patches is hereafter referred to as a "ring-around." The 
target will have a group of clusters of patches, including the neutral patch and the 
deviated patches. Each cluster is used to produce one triad of red, green and blue 
drive code values that, Vhen used together produce one neutral density patch. The 
ultimate objective is that\ach group of clustered patches is used to find a single 
triad of red, green and blue chive code values that produce a single aim neutral. 
Then combining all such solutions for all the clusters of patches on the target 
perform separate curvefits on theVdensity vs. code value solutions for each of the 
color channels as if they were independent: 

D r = f rcd (dCV r )\ 

D g = WdCV^X 

D b = f b iue(dCV b ), \ 
with the understanding that these equations arte valid only when used together and 
the collective solutions produce a triad of densities that correspond to an aim 
neutral. This set of functions can be called the neuftal characteristic curves of the 
rendering device because they describe its native behavior when producing neutral 
colors. \ 

These 3 curvefit functions are then combined with the three aim 
curves to produce new calibration look up tables as follows. The aim curves map 
neutral input code values to sets of densities producing neutral. The neutral 
characteristic curves map the drive code values to densities producing neutral. 
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Thus a new mapping between input code values and output code values can be 
determined using the neutral densities as the common link. It is to be understood 
that all of the steps described above are preferably performed automatically in a 
digital image processor, with the possible exception of manually transferring 
rendered test patches from the reproduction apparatus to a scanner. 

According to one preferred embodiment of the present invention, a 
set of, say, one to two dozen neutral points are selected along the neutral axis (i.e., 
the locus of neutral sets-of-densities) in a three-dimensional density space whose 
axes corresponding to the red, green and blue densities. For each neutral point, 
four non-neutral points are selected such that the four non-neutral points sit on a 
common plane wherein (1) the neutral axis is normal to the common plane, (2) the 
common plane is further out (darker) on the neutral axis than the corresponding 
neutral point, (3) the four non-neutral points represent the corners of a square in 
the common plane such that the square is centered on the point where the neutral 
axis cuts the common plane, and (4) all four non-neutral points are roughly 
equidistant from each other. 

The invention has been described in detail with particular reference 
to preferred embodiments thereof, but it will be understood that variations and 
modifications can be effected within the spirit and scope of the invention. 



